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Abstract 



In Randall-Sundrum-type brane-world cosmologies, the dynamical equa- 
tions on the three-brane differ from the general relativity equations by terms 
that carry the effects of imbedding and of the free gravitational field in the five- 
dimensional bulk. Instead of starting from an ansatz for the metric, we derive 
the covariant nonlinear dynamical equations for the gravitational and matter 
fields on the brane, and then linearize to find the perturbation equations on 
the brane. The local energy-momentum corrections are significant only at 
very high energies. The imprint on the brane of the nonlocal gravitational 
field in the bulk is more subtle, and we provide a careful decomposition of this 
effect into nonlocal energy density, flux and anisotropic stress. The nonlocal 
energy density determines the tidal acceleration in the off-brane direction, 
and can oppose singularity formation via the generalized Raychaudhuri equa- 
tion. Unlike the nonlocal energy density and flux, the nonlocal anisotropic 
stress is not determined by an evolution equation on the brane, reflecting 
the fact that brane observers cannot in general make predictions from initial 
data. In particular, isotropy of the cosmic microwave background may no 
longer guarantee a Friedmann geometry. Adiabatic density perturbations are 
coupled to perturbations in the nonlocal bulk field, and in general the system 
is not closed on the brane. But on super-Hubble scales, density perturba- 
tions satisfy a decoupled third-order equation, and can be evaluated by brane 
observers. Tensor perturbations on the brane are suppressed by local bulk 
effects during infiation, while the nonlocal effects can serve as a source or a 
sink. Vorticity on the brane decays as in general relativity, but nonlocal bulk 
effects can source the gravito-magnetic field, so that vector perturbations can 
be generated in the absence of vorticity. 
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I. INTRODUCTION 



Einstein's theory of general relativity breaks down at high enough energies, and is likely 
to be a limit of a more general theory. Recent developments in string theory indicate 
that gravity may be a truly higher- dimensional theory, becoming effectively 4-dimensional 
at lower energies. These exciting theoretical developments need to be accompanied by 
efforts to test such higher-dimensional theories against their cosmological and astrophysical 
implications. In that spirit, we investigate here a particular class of models, showing how 
their dynamical properties generalize those of Einstein's theory, and discussing the broad 
implications of these generalizations for cosmological dynamics. 

In many higher- dimensional gravity theories inspired by string theory, the matter fields 
are confined to a 3-brane in 1 -|- 3 -|- o? dimensions, while the gravitational field can propagate 
also in the d extra dimensions (see, e.g., |[^). It is not necessary for the d extra dimensions 
to be small, or even finite: recently Randall and Sundrum |]^ have shown that for d = 1, 
gravity can be localized on a single 3-brane even when the fifth dimension is infinite (see 
also 0). An elegant geometric formulation and generalization of the Randall- Sundrum 
scenario has been given by Shiromizu, Maeda and Sasaki The Friedmann equation on 
the brane in these models is modified by both high-energy matter terms and a term carrying 
nonlocal bulk effects onto the brane. The Friedmann brane models have been extensively 
investigated (see, e.g., 0-0]), and inflationary scalar perturbations in these models have also 
been considered [|]. The models are compatible with observations subject to reasonable 
constraints on the parameters. A broader study of cosmological dynamics, i.e., for induced 
metrics more general than the simple Friedmann case, has not been done. In particular, the 
analysis of perturbed Friedmann models also remains to be done. (Considerable work has 
been done on perturbations of flat brane metrics; see, e.g., p|j9|JT0|l.) 

In this paper, we initiate a study of nonlinear and perturbed cosmological dynamics in 
Randall-Sundrum-type brane-world models, generalizing some important results in general 
relativity. We find the bulk corrections to the propagation and constraint equations, us- 



ing the covariant Lagrangian approach ITjjIil]- This approach is well-suited to identifying 
the geometric and physical quantities that determine inhomogeneity and anisotropy on the 
brane, and it is also the basis for a gauge-invariant perturbation theory Our first task 



is to identify and interpret the covariant physical content of the bulk effects on the brane. 
Local effects lead to quadratic corrections of the density, pressure and energy flux. The 
nonlocal effects of the free gravitational field in the bulk are transmitted by a Weyl projec- 
tion term, which we decompose into energy density, energy flux and anisotropic stress parts. 
We calculate the gravitational (tidal) and non-gravitational acceleration of fluid world-lines, 
finding the role of the nonlocal energy density in localization of gravity, and showing how 
the world-lines have a non-gravitational acceleration off the brane at high energies. During 
inflation, the acceleration is directed towards the brane. 

We derive the propagation ('conservation') equations governing the nonlocal energy den- 
sity and flux parts; the evolution of the anisotropic stress part is not determined on the 
brane. These nonlocal terms also enter into crucial dynamical equations, such as the Ray- 
chaudhuri equation and the shear propagation equation, and can lead to important changes 
from the general relativistic case. For example, it is possible via the nonlocal term to avoid 
the initial singularity in a nonrotating model without cosmological constant. Nonlocal ef- 
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fects also mean that isotropy of the cosmic microwave background (CMB) may no longer 
guarantee a Friedmann geometry. 

The covariant nonlinear equations lead to a covariant and gauge-invariant description of 
perturbations on the brane. We derive the equations governing adiabatic scalar perturba- 
tions, which are not in general closed on the brane, because of nonlocal effects. However, 
on super-Hubble scales, the density perturbations satisfy a decoupled third-order equation, 
with an additional nonlocal degree of freedom, and can therefore be evaluated by brane ob- 
servers. Tensor perturbations cannot be determined by brane observers on any scales. The 
local bulk effects tend to enhance tensor perturbations during non- inflationary expansion 
and suppress them during inflation. Nonlocal bulk effects can in principle act either way. 
Vorticity on the brane decays as in general relativity, but bulk effects act as a source for the 
gravito-magnetic field, and hence vector perturbations, on the brane. 

Our results remain incomplete in one fundamental aspect, i.e., we do not provide a 
description of the gravitational field in the bulk, but confine our investigations to effects 
that can be measured by brane-observers. In order to fill this gap, one would need to study 
the off-brane derivatives of the curvature, which are given in general in This is an 

important topic for further research. 

The 5-dimensional field equations are Einstein's equations, with a (negative) bulk cos- 
mological constant A and brane energy-momentum as source: 



Gab — 



-kgAB + 5{x){-^9AB + TAB} ■ (1) 



The tildes denote the bulk (5-dimensional) generalization of standard general relativity quan- 
tities, and = Svr/Mp, where Mp is the fundamental 5-dimensional Planck mass, which 
is typically much less than the effective Planck mass on the brane, Mp = 1.2 x 10^^ GeV. 
The brane is given by x = 0, so that a natural choice of coordinates is = (x^,x), 
where x'^ = {t, x^) are spacetime coordinates on the brane. The brane tension is A, and 
9ab = dAB — nATiB is the induced metric on the brane, with ua the spacelike unit normal to 
the brane. Matter fields confined to the brane make up the brane energy-momentum tensor 
Tab (with TABn^ = 0). 

Although it is usually assumed that the spacetime is exactly anti-de Sitter in the absence 
of a brane (A = = Tab), this is not necessarily the case. The brane-free bulk metric can 
be any solution of Gab = —^'^^Qabi including non-conformally fiat solutions. When the 
induced metric on the brane is Friedmann, then the 5-dimensional Schwarzschild-anti-de 
Sitter metric is a solution of Eq. ([I|) |^ . However, more general bulk metrics are in principle 
possible. 

The field equations induced on the brane are derived via an elegant geometric approach 
by Shiromizu et al. 0], using the Gauss-Codazzi equations, matching conditions and Z2 
symmetry. The result is a modification of the standard Einstein equations, with the new 
terms carrying bulk effects onto the brane:[| 

G ^1/ = A(7^jy + K T^y + K S^i, S^i, , (2) 



^ For clarity and consistency, we have changed the notation of 
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where = Sir/M^. The energy scales are related to each other via 



A = 6— , A = ^ 



3Mp3 



(3) 



The bulk corrections to the Einstein equations on the brane are of two forms: firstly, 
the matter fields contribute local quadratic energy-momentum corrections via the tensor 
S^u, and secondly, there are nonlocal effects from the free gravitational field in the bulk, 
transmitted via the projection S^i, of the bulk Weyl tensor. The matter corrections are given 
by 

c J_T '^T It T" -I Ln I'^T T<^P (T ""i^ (A\ 

The projection of the bulk Weyl tensor is0 

£ab = CACBon^n^ , (5) 

which is symmetric^ and traceless {£[ab] = = £a^) and without components orthogonal 
to the brane, so that Sab^^ = and Sab S^,^6a'^Sb'^ as x ^ 0. The Weyl tensor Cabcd 
represents the free, nonlocal gravitational field in the bulk. The local part of the bulk 
gravitational field is the Einstein tensor Gab, which is determined via the bulk field equations 
(0). Thus S^u transmits nonlocal gravitational degrees of freedom from the bulk to the brane, 
including tidal (or Coulomb), gravito-magnetic and transverse traceless (gravitational wave) 
effects. 

II. COVARIANT DECOMPOSITION OF BULK EFFECTS 

We now provide a covariant decomposition of the bulk correction tensors given by Shi- 
romizu et al. 0]. 

A. Local bulk effects 

For any matter fields (scalar fields, perfect fiuids, kinetic gases, dissipative fiuids, etc.), 
including a combination of different fields, the general form of the brane energy-momentum 
tensor can be covariantly given as 

Tf,^ = pUfj,Uy + phf,^ + vr^^ + 2q^^,Uy) . (6) 



^ This projection is called the 'electric' part of the Weyl tensor in Q, but the term is potentially 
misleading, since the electric part is associated with projection on a timelike vector |11,14|, and 



n is spacelike. should not be confused with the electric part of the brane Weyl tensor, Efj^^, 
defined below. 

3 (Square) round brackets enclosing indices denote (anti-)symmetrization. 
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The decomposition is irreducible for any chosen 4-velocity u^. Here p and p are the energy 
density and isotropic pressure, and h^^ = g^^ + u^u^ projects orthogonal to u'^. The energy 
flux obeys = g^^), and the anisotropic stress obeys n^^ = ^{^u), where angle brackets 
denote the projected, symmetric and tracefree part: 



Equations (H) and (W) imply the irreducible decomposition 



S, 



J_ 

24 



h 



flU 



For a perfect fluid or minimally-coupled scalar field 

S„ 



(7) 



in agreement with [^]. The quadratic energy- momentum corrections to standard general 
relativity will be significant for H'^p^ ^ K^p, i.e., in the high-energy regime 



B. Nonlocal bulk effects 

The symmetry properties of £^j,i, imply that in general we can decompose it irreducibly 
with respect to a chosen 4-velocity field as 



S 



U (u^u, + I V) + ^/^^ + 2Q(mM.)] • (8) 



The factor (k/k)^ is introduced for dimensional reasons. Here 

U = - [^^^^.u^u" 

is an effective nonlocal energy density on the brane, arising from the free gravitational field 
in the bulk. This nonlocal energy density need not be positive (see below). There is an 
effective nonlocal anisotropic stress 

on the brane, arising from the free gravitational field in the bulk, and 

4 



Q 



is an effective nonlocal energy flux on the brane, arising from the free gravitational field in 
the bulk. 



5 



If the induced metric on the brane is flat, and the bulk is anti-de Sitter, as in the original 
Randall-Sundrum scenario 0], then 8^,^ = 0. Treating this as a background, it follows that 
tensor (transverse traceless) perturbations on the brane arising from nonlocal bulk degrees 
of freedom are given by 

U = 0=Q^, B'V^, = , (9) 
where is the totally projected part of the brane covariant derivative: 

D^F°-...^ = V^"7---VV.F^-...5. 

Equation (^) provides a covariant characterization of brane tensor perturbations on an anti- 
de Sitter background. 

In cosmology, the background induced metric is not flat, but a spatially homogeneous 
and isotropic Friedmann model, for which 

I)^U = Q^ = V^, = 0. (10) 

Thus for a perturbed Friedmann model, the nonlocal bulk effects are covariantly and gauge- 
invariantly described by the flrst-order quantities D^U, Q^, Vfj,u- Since W 7^ in general in 
the Friedmann background |^,|^, it follows that for tensor perturbations on the brane 

B^U = 0=Q^, D''V^, = 0. (11) 

Scalar perturbations (Coulomb-like bulk effects) will be characterized by 

Qm = D;.Q, P;.. = D(^D,)P, (12) 

for some scalars Q and V, while for vector perturbations (gravito- magnetic-like bulk effects) 

D^Q^ = , V^, = D(^P,) , D^P^ = . (13) 



C. Gravitational and non-gravitational acceleration 

In order to flnd the role of bulk effects in tidal acceleration on the brane, we start from 
the relation 

SavU^v!" = limC-ABCDU'^n^u^n^ , 

where is an extension off the brane of the 4- velocity (with u'^riA = 0). The tidal accel- 
eration in the direction measured by comoving observers is —fiAR^BCDU^n'^u^- Now 

RaBCD = CaBCD + I {gA[cRD]B + ^B[Z)-Rc]a} " \RgA[cgD]B , 

so that by the field equation (|l]) (and recalling that Tab^^ = 0), 

-RABCon^u^rfu^ = -Sabu^u^ + \k^K . 
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Thus the comoving tidal acceleration on the hrane, in the off-brane direction, is 



K 



K 



4 



- U + lU'A. (14) 



Since A is negative, it contributes to acceleration towards the brane. This reflects the 
confining role of the negative bulk cosmological constant on the gravitational field in the 



warped metric models of Randall- Sundrum type. Equation (|T^ shows that localization of 
the gravitational field near the brane is enhanced by negative U, while positive U acts against 
confinement, by contributing to tidal acceleration away from the brane. This picture is 
consistent with a Newtonian interpretation, where the gravitational field carries negative 
energy density. 

The non-gravitational acceleration of fiuid world-lines on the brane is 

= u^Vbu^ at X = . 
Locally, near the brane, the metric may be written as 

ds'^ = dx^ + g^iuix"", x)dx^dx^ , 

so that 

fV(a;,0) = rV(a;)5a^ - |^?^.,x(x,0)5/. 

This allows us to express the 5-dimensional acceleration in terms of the 4-dimensional 
acceleration A^ = u^SJ yU^ on the brane. First, we use the covariant form g^v.x — l^-a.9ii.v, 
where is the Lie derivative along n^. Then we use the expression for the extrinsic 
curvature of the brane [01 



which leads to 

A^{x, 0+) = A''{x)6/ - K+^{xy\x)u''{x)n'^{x, 0+) , (15) 

on the brane. The extrinsic curvature is given in terms of the brane tension and energy- 
momentum by 

= -l^"^ [^9i,u + ^T^y + (p - 3p)^^^] . 



Substituting in Eq. (p!5D, we find 



A^{x, 0+) = A>'{x)5^^ + [2p{x) + 3p(x) - A] n^(x, 0+) . (16) 

It follows that tlaA^ is nonzero on the brane, i.e., there is a non- gravitational acceleration 
of fluid world-lines orthogonal to the brane. The direction depends on the sign of 2p+3p — A : 
for2p-\-3p — X > 0, the transverse acceleration is away from the brane, while for 2p-\-3p— X < 
0, it is towards the brane. If the pressure is positive, then at high energies 2p + 3p — X > 0, 
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so that either the brane must accelerate, or there must be a non-gravitational mechanism 
for keeping matter on the brane. 

Inflationary expansion, when pressure is negative, can change this situation. Inflation 
on the brane is characterized by 



f<- (17) 




and this condition implies that 2p + 3p — X < 0. Thus during inflation on the brane, the 
transverse acceleration is towards the brane; inflation acts as a non- gravitational mechanism 
keeping matter on the brane. 



D. Effective total energy-momentum tensor 

All the bulk corrections may be consolidated into effective total energy density, pressure, 
anisotropic stress and energy flux, as follows. The modified Einstein equations take the 
standard Einstein form with a redefined energy-momentum tensor: 



where 



Then 



G 



T. 



Tun -\- „ S nu „ £u 



tot 



P 



tot 



P+ — 



24 

.4 



IT 



Jot 



K 



,tot 



(18) 
(19) 

(20) 
(21) 
(22) 
(23) 



(Note that / is dimensionless.) 

These general expressions simplify in the case of a perfect fluid (or minimally coupled 
scalar field, or isotropic one-particle distribution function), i.e., for = = vr^j,. However, 
the total energy flux and anisotropic stress do not vanish in this case in general: 



K 



9 



TT 



tot 



K 



-O 

5 ^^A* ' ^iv 



K 



V 



Thus nonlocal bulk effects can contribute to effective imperfect fluid terms even when the 
matter on the brane has perfect fluid form. 
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III. LOCAL AND NONLOCAL CONSERVATION EQUATIONS 

As a consequence of the form of the bulk energy-momentum tensor in Eq. (|I|) and of 
Z2 symmetry, it follows that the brane energy-momentum tensor separately satisfies the 
conservation equations, i.e., 

V%. = . (24) 

Then the Bianchi identities on the brane imply that the projected Weyl tensor obeys the 
constraint 

V^^^. = K^V^S^, , (25) 

which shows that its longitudinal part is sourced by quadratic energy-momentum terms, 
including spatial gradients and time derivatives. Thus evolution and inhomogeneity in the 
matter fields generates nonlocal Coulomb-like gravitational effects in the bulk, which 'back- 
react' on the brane. 

The brane energy-momentum tensor and the consolidated effective energy-momentum 
tensor are both conserved separately. Conservation of T^y gives the standard general rela- 



tivity energy and momentum conservation equations 15 



p + e (p + p) + D^g^ + iN'q^ + a'^^TT^, = , (26) 
+ 1©?/. + D^p + (p + v)K + D'^TT^^ + A'^vr^^ + o^X - g]/. = . (27) 

A dot denotes m'^Vi,, = D'^u^ is the volume expansion rate of the congruence, = 
= ^(/i> is its 4-acceleration, a^,^ = Dd^Uu) is its shear rate, and cu^ = — |curlM^ = cj^^) is 
its vorticity rate. The covariant spatial curl is given by |T2| 

curl = e^^p-D'^V^ , curl W,, = £«^(^D"W^^) , 

where e^^o- is the projected alternating tensor. The covariant cross-product is 

The conservation of T*°* gives, upon using Eqs. (^OD-(^), a propagation equation for 
the nonlocal energy density lA: 

U + IQU + D'^Q^ + 2A^Q^ + a^'^P^, 

= [67i^''7r^, + 6(p + p)a^''7c^, + 26 i2q^q^ - ti^^ti^,) 
+ 2A^'q''n^, + 4g^D^p + q^W-n^, + Tr^^D^g, - 2a^'^7r„^7r,° - 2a^'''q^q,] , (28) 

and a propagation equation for the nonlocal energy flux Q^: 

Q(M> + feS^ + \D^U + pA^ + D^V^, + A'V^, + a^.Q" - [uj, Q]^ 



24'^ 



-4(p + p)D^p + g^TT^^,) + 7r/D,(2p + hp) + 6(p + p)D^ 



- e^^pu^TiC^q, + 4(p + p)Qq^ + Qq^A'q, + f A^g^g. + 4g,a°^7r„ J . (29) 
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These equations may be thought of as conservation equations on the brane for the nonlocal 
energy density and energy flux due to the free gravitational field in the bulk. In general, 
the 4 independent equations determine 4 of the 9 independent components of E^^ on the 
brane. What is missing, is an evolution equation for V^y (which has up to 5 independent 
components). Thus in general, the projection of the 5-dimensional field equations onto the 
brane does not lead to a closed system. Nor could we expect this to be the case, since 
there are bulk degrees of freedom whose impact on the brane cannot be predicted by brane 
observers. These degrees of freedom could arise from propagating gravity waves in the bulk, 
possibly generated by inhomogeneity on the brane itself. The point is that waves which 
penetrate the 5th dimension are governed by off-brane bulk dynamical equations. Our 
decomposition of E^y has shown that the evolution of the nonlocal energy density and flux 
(associated with the scalar and vector parts of £^y) is determined on the brane, while the 
evolution of the nonlocal anisotropic stress ( associated with the tensor part of £^u) is not. 
If the nonlocal anisotropic stress contribution from the bulk field vanishes, i.e., if 

P^. = , 

then the evolution of E^y is determined by Eqs. ( ]28|) and (pOf). A special case of this arises 
when the induced metric on the brane is Friedmann, i.e., when Eq. ([ToD holds. Then 
has only 1 independent component, W, and it is determined by Eq. (^) (see below), with 
Eq. (^9]) reducing to = 0. Another case when the equations close on the brane is when 
the brane is static (see |jIB|). 

In the perfect fluid case, the conservation equations (^)- (|29|) reduce to: 

p + e(p + p) = o, (30) 

'D^.P+{p + p)A^ = Q. (31) 

For a minimally-coupled scalar fleld, p = |V/_iV^V'^v' + V{ip) and p = ^V^ipV^ip — V{{p). In 
the adiabatic case, Eq. (|3TD gives 

A, = -^D,p, cl = ^. (32) 
p + p p 



The nonlocal conservation equations (^) and (|29D reduce to 



W + few + D^Q,, + 2A^Q,, + a^^'V^, = , (33) 
k\p + p)D,p. (34) 



1-4 
"6 



Equation (0) shows that if Sfj_i, = and the brane energy-momentum tensor has perfect 
fluid form, then the density p must be homogeneous The converse does not hold, i.e., 
homogeneous density does not in general imply vanishing S^^. This is readily apparent from 
Eq. (|34D . A simple example is provided by the Friedmann case: Eq. (|3^) is trivially satisfled, 
while Eq. (^) gives the 'dark radiation' solution 

«^«„(^)\ (35) 
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A simple generalization of the Friedmann case is the purely Coulomb-like case, = 



V^iTj. Equation (35) still holds, but with a an average scale factor, which is in general 
inhomogeneous. Equation (0) reduces to a constraint on the acceleration. Local momentum 
conservation already provides the constraint in Eq. (p2D. It follows that in the purely 
Coulomb-like case, the spatial gradient of the nonlocal energy density must be proportional 
to that of the local energy density: 



clu - K^{p + py 

2(p + p) 



Linearization about a Friedmann background does not change Eqs. (|30|) and (0), but 
Eqs. (PD and (|3l) lead to 



W + |eW + D'^Q^ = 0, (36) 
+ 4HQ^ + iD^W + pA^ + D'V^, = -|/€^(p + p)D^p , (37) 

where H is the Hubble rate in the background. The nonlocal tensor mode, which satisfies 
D^V^y = 07^ V^y, does not enter the nonlocal conservation equations. 

IV. PROPAGATION AND CONSTRAINT EQUATIONS 



Equations (p6D-(p^) are propagation equations for the local and nonlocal energy density 
and flux. The remaining covariant equations on the brane are the propagation and constraint 
equations for the kinematic quantities and the free gravitational field on the brane. The 
kinematic quantities govern the relative motion of neighboring fundamental world-lines, and 
describe the universal expansion and its local anisotropies. The locally free gravitational 
field on the brane is given by the brane Weyl tensor C ^1,^/3 ■ This splits irreducibly for a given 

into the gravito-electric and gravito-magnetic fields on the brane: 



where E^y must not be confused with £^y. The Ricci identity for and the Bianchi identities 
= V [^{—Ru]a + \R9u]a) producc the fundamental evolution and constraint equa- 
tions governing the above covariant quantities Einstein's equations are incorporated via 



the algebraic replacement of the Ricci tensor R^^ by the effective total energy-momentum 
tensor, according to Eq. (pISl). These are derived directly from the standard general relativ- 
ity versions (see, e.g., [I3|) by simply replacing the energy-momentum tensor terms p, ... by 



p , .... The result for a general imperfect fluid is given in Appendix A. 

A. Nonlinear equations 

For a perfect fluid or minimally-coupled scalar field, the equations in Appendix A reduce 
to the following. 
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Expansion propagation (generalized Raychaudhuri equation): 



1 



12 \k 
Vorticity propagation: 



n^p{2p + 3p) + 12W 



1 fK 



Shear propagation: 

Z \ AC 

Gravito-electric propagation: 



E{f,^) + e^^i, - curl H^i, + {p + p)a^ 

Gravito-magnetic propagation: 

H{,^u) + + curl^^^i, — Saai^H^)"" + u'^EafsifiHiy/ + 2A°'eap(^E^)^ 

4 



1 



~ 2 

Vorticity constraint: 
Shear constraint: 



curiP^^ - ?>uj(^Qu) + (7°' {^e^)apQ!' 



DV^, - curla;^ - fD^O + 2[a;, A\^ = - j . 

Gravito-magnetic constraint: 

curl (7^^ + D(^a;j.) - if^^ + 2Ai^^uj^) = . 
Gravito-electric divergence: 

D'^E^, - Ik^D^p - [a, HI + SH^^cu" 
- ^ (^) {«^VD/.P + 6D^W - 692^ - QD'-T'^, + Qa^.Q'^ - 27[^, Q]^ . 
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Gravito-magnetic divergence: 



- (-) {k*p{p + p)uj^ - 3curl + 8Wcu^ - 3[cr, V\ - SV^.u"} . (47) 



6 \K, 

Here the covariant tensor commutator is 

The standard 4-dimensional general relativity results are regained by setting all right hand 
sides to zero in Eqs. (|38D-(^7D. 



Together with Eqs. (PD|)-(^), these equations govern the dynamics of the matter and 
gravitational fields on the brane, incorporating both the local (quadratic energy-momentum) 
and nonlocal (projected Weyl) effects from the bulk. These effects give rise to important 
new driving and source terms in the propagation and constraint equations. The vorticity 
propagation and constraint, and the gravito-magnetic constraint have no explicit bulk effects, 
but all other equations do. Local and nonlocal energy density are driving terms in the 
expansion propagation, and note that these are the terms that determine the gravitational 
and non-gravitational acceleration transverse to the brane. The spatial gradients of local 
and nonlocal energy density provide sources for the gravito-electric field. The nonlocal 
anisotropic stress is a driving term in the propagation of shear and the gravito-electric/ 
-magnetic fields, and the nonlocal energy fiux is a source for shear and the gravito-magnetic 
field. 

In general the system of equations is not closed: there is no evolution equation for the 
nonlocal anisotropic stress V^u, which carries the tensor modes in perturbed solutions. If 
we set Sij,iy = to close the system, i.e., if we allow only local matter effects from the fifth 
dimension, then, as noted above, the density is forced to be homogeneous. This is clearly 
too restrictive. A less restrictive way of closing the system is to assume that the nonlocal 
anisotropic stress vanishes, i.e., V^u = 0. However, this rules out tensor modes arising from 
the free field in the bulk, and also limits the scalar and vector modes, which will in general 
also be present in Vfj.u- 



B. Linearized equations 

We have derived the exact nonlinear equations that govern gravitational dynamics on 
the brane as seen by brane observers. These equations hold for any geometry of the brane, 
and they are fully covariant on the brane. In particular, this means that we can linearize the 
equations by taking a suitable limit, and without starting from a given background solution. 
In this way we avoid the need for choosing coordinates, and we deal directly with covariant 
physical and geometric quantities, rather than metric components. 

The limiting case of the background Friedmann brane is characterized by the vanishing 
of all inhomogeneous and anisotropic quantities. These quantities are then first-order of 
smallness in the linearization scheme, and since they vanish in the background, they are 



gauge-invariant |13|. The standard general relativity scheme is modified by the additional 
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degrees of freedom arising from bulk effects. In particular, the generalized Friedmann equa- 
tion on the brane is UBI 



K 



1 



3 \K 



(48) 



where = 0, ±1. 

Linearization about a Friedmann brane model of the propagation and constraint equa- 
tions leads to the reduced system: 



+ ie' - DM^ + W{p + 3p) - A 



1 fH 



12 \k 



K^p{2p + 3j9) + 12U 



LUfj, + 2Huj^ + icurl = , 



1 



2 V/t 



E^^ + 3iJE^i, - curl H^^ + \k^{p + p)a^y 



1 /K 



12 

if^,, + ?>HH^y + curl E^,, 
D'^u;^ = , 

curl a^y + D(^u;^) - ff^^ 

1 ,.2i 



\ (- ] curlP, 



2 



-) 



4 

) Q 





0, 






'K 


r 






1 




6 





{kV(p + p) + 8U} uj^ - 3curl 



(49) 
(50) 

(51) 

(52) 

(53) 
(54) 

(55) 
(56) 

(57) 
(58) 



These equations, together with the linearized conservation equations (^), (|3ll) , ( ^6|) and 
(^), are the basis for a covariant and gauge-invariant description of perturbations on the 
brane. The local bulk effects (quadratic energy-momentum effects) are purely scalar, as 
is the nonlocal energy density. The nonlocal energy flux has in general scalar and vector 
modes: 



and scalar, vector and tensor modes enter the nonlocal anisotropic stress: 



(59) 



(60) 



In these equations, an overbar denotes a transverse (divergence-free) quantity: 



0. 



14 



V. NONLINEAR AND PERTURBATIVE DYNAMICS 



Bulk effects introduce new degrees of freedom into the dynamics on the brane, subject to 
the additional nonlocal 'conservation' equations ( |33D and (Q). Standard results in general 
relativity may or may not continue to hold under this higher-dimensional modification. We 
now use the conservation, propagation and constraint equations to generalize some standard 
results of 4-dimensional general relativity. 



A. CMB isotropy and brane homogeneity 



In standard general relativity, the isotropy of the CMB radiation has crucial implications 
for the homogeneity of the universe. If all fundamental observers after last scattering observe 



an isotropic CMB, then it follows from a theorem of Ehlers, Geren and Sachs that 
the universe must have a homogeneous Friedmann geometry. This has been generalized 
to the almost-isotropic case |T^, providing a foundation for the perturbative analysis of 
CMB anisotropics. The Ehlers-Geren-Sachs theorem is based on the coUisionless Boltzmann 
equation, and on the kinematic/dynamic characterization: 



A, 



o-^i. = 7^ and 







Friedmann geometry. 



(61) 



Bulk effects do not change the Boltzmann equation, but they do mean that the Friedmann 
characterization is no longer in general true on the brane. While the gravito-magnetic 



constraint Eq. 



still leads to Hn^, = 0, the shear propagation equation 



no longer 



forces Efj^i, = 0, because of the nonlocal term V^u, so that the intrinsic metric need not be 
conformally flat. 

A consistent solution on the brane of the system of (nonlinear) conservation, propagation 
and constraint equations can be given as follows. We take 

D^p = D^p, = D,,W = D^0 = O, Q^ = = D^P^,. 

Then the system of equations reduces to the consistent set 



P+Qp 

Pr + |0Pr 

w + few 



0, 
0, 
0, 



e + + y{p + 2pr) - A = ( - ) \^\p + p,)(2p + 3p,) + uu 



V, 



In general, provided that the propagation equation for V^v is consistent with the bulk 
geometry, E^y need not be zero, so that the brane geometry need not be Friedmann, although 
it is asymptotically Friedmann. 

Thus it is in principle possible via nonlocal bulk effects that isotropic CMB does not 
force the brane metric to be Friedmann. 
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B. Generalized gravitational collapse 



Another important question is how the higher-dimensional bulk effects modify the picture 
of gravitational collapse/ singularities, which depends on Raychaudhuri's equation ||lT| . 



The generalized Raychaudhuri equation (^) governs gravitational collapse and initial 
singularity behavior on the brane. The local energy density and pressure corrections, 

^^V(2p + 3p) , 

further enhance the tendency to collapse, if 2p + 3p > 0. This condition will be satisfied 
in thermal collapse (or time-reversed expansion), but it is violated during very high-energy 
inflation (p ^ A), by Eq. (|T^, and in that case the local bulk term acts to further accelerate 
expansion. This is consistent with the results given in ^j. Thus local bulk effects at high 
energy reinforce the formation of singularities during thermal collapse, as predicted in general 
relativity, while further accelerating expansion during high- energy inflation. 
The nonlocal term. 



K 



can act either way, depending on its sign. As shown from Eq. (p!^, a negative U enhances 
the localization of the gravitational field on the brane. In this case, the effect of U is to 
counteract gravitational collapse. A positive U acts against localization, and also reinforces 
the tendency to collapse. 

If higher-dimensional corrections to Einstein's theory tend to prevent singularities, then 
the effective energy density U on the brane of the free gravitational field in the bulk should 
be negative. In this case, U also acts to reinforce confinement of the gravitational field to 
the brane. 



C. Cosmological scalar perturbations 

The linearized equations on the brane derived in the previous two sections encompass 
scalar, vector and tensor modes. In order to covariantly (and locally) separate out the 
scalar modes, we impose the condition that all perturbative quantities are spatial gradients 
of scalars, i.e.. 

The identities in Appendix B, the vorticity constraint equation ( p^ and the gravito-magnetic 
constraint equation (|56|) then show that 

curl = = curl W^, , B^W^, = p^iB^W) , uj^ = = H^,, (62) 

as in standard general relativity. 

If we choose the fundamental 4-velocity such that cr^jy = 0, which is the covariant 
analogue of the longitudinal or conformal Newtonian gauge in metric-based perturbation 
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theory (see |lT9| , pO| for further discussion), then the shear propagation equation (|51|) becomes 
a constraint determining the brane tidal tensor: 



1 /K 



2 \k 



V 



Here $ is the relativistic generahzation of the gravitational potential, defined by = D^<l>, 
and V is the potential for the nonlocal anisotropic stress, defined in Eq. (160|) . It follows that 
nonlocal bulk effects lead to a change in the gravitational tidal potential (in longitudinal-like 
gauge): 



1 /K 



2 \k 



V. 



(63) 



In the general case, i.e., when is not chosen to give a^^ = 0, this simple relation does not 
hold. 

In order to derive the equations governing density perturbations in the general case, we 
define the density and expansion gradients (as in 



A 



P 



^ - -D^p , = aD„6 



(64) 



and the (dimensionless) gradients describing inhomogeneity in the nonlocal quantities: 

1 



p p 



ap 



(65) 



where Q is defined in Eq. (0). Then we take the spatial gradient of the energy conservation 
equations (|30|) and (|36D and the generalized Raychaudhuri equation (^9|), using the identities 
in Appendix B and the adiabatic equation (^2]). We arrive at the following system of 
equations: 



A„ = 3wHA, 



-2HZ„ 



l + w 



-) pu. 



K 



+ |k^(4 + 3w)p - 



% = {3w - 1)HU, 



4c^ 
l + w 



.p. 



K 

HA, 



2d 



U 



= (1 - 3w)EQ^ - —U^ - laD'P^ + — 



1 

6a 



1 + w J p 

-S) 



Z^ - aD^Q^ , 



--K^p\l+w) 



A 



(66) 

(67) 
(68) 
(69) 



1 + w J p 

where w = p/p. The background Friedmann equation (^) relates H and p, with U in the 
background given by Eq. (|35D . 

In standard general relativity, only the first two equations apply, with k, set to zero in 
Eq. (|67D . In this case we can decouple the density perturbations via a second-order equation 
for A^, whose independent solutions are adiabatic growing and decaying modes. Nonlocal 
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bulk effects introduce important changes to this simple picture. First, we note that there is 
no equation for P^, so that in general, scalar perturbations on the brane cannot be predicted 
by brane observers without additional information from the unobservable bulk. 

However, there is a very important exception to this, arising from the fact that and 
only occur in Eqs. (|66|)-(p8|) via the Laplacian terms D^Q^ and D^P^, and the latter 
term is the only occurrence of in the system. It follows that on super-Hubble scales, the 
system does close, and brane observers can predict scalar perturbations from initial conditions 
intrinsic to the brane. The system reduces to 3 coupled equations in A^, and [/^, plus 
an equation for Q^, which is determined once the other 3 quantities are solved for. One can 
decouple the density perturbations via a third-order equation for A^. The nonlocal energy 
density plays the role of a non-interacting radiation fluid with the same velocity as the 
ordinary fluid, and inhomogeneity in U introduces an additional entropy-like scalar mode. 
As may have been expected, this additional mode is absent during radiation-domination; in 
this case Eqs. (pO) and (pSI) show that 



U„ 



Po 



In principle, it is straightforward to solve the coupled system on super-Hubble scales, 
although numerical techniques will be necessary. We do not attempt particular solutions 
here. However, it may be instructive to see the decoupled third-order equation for density 
perturbations during matter-domination, on a fiat background: 



A, 



2HA^ + 
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7 2 2~4 2 



3^^ 



A^ + 



2~4 
3K p 



pHA^ = . (70) 



D. Cosmological tensor perturbations 

Tensor perturbations are covariantly characterized by 



0, Ti^'W, 







where / = p,p,Q,U, and W^i, = a^y, E^y, H^^^V^i,. Then all the conservation equations 
reduce to background equations, and the system of linearized propagation and constraint 
equations in the previous section reduces to: 







E;, + she;, - cm\H^, + \k^{p + p)a^. 



K 



12 \k 

H^, + 3HH^, + cnTlE: 



{/€V(P + P) + 8W} c^M- + 12 {V,.u + 2HP, 
0, 



curler 



H, 







(71) 



(72) 

(73) 
(74) 



where the overbars denote transverse tensors, and 
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Since there is no equation for V^u-, the system of equations does not close on the brane: brane 
observers cannot evaluate tensor perturbations on the brane without additional information 
from the unobservable bulk. This remains true on super-Hubble scales, unlike the scalar 
perturbation case. 



Equations (^) and ( |74|) show that the shear is a gravito-potential for E* and 



Using the identities in Appendix B, we can derive the following covariant wave equation for 
the shear: 



4 



2A + iK'(p-3j9) 



cr 



(75) 



For adiabatic tensor perturbations in standard general relativity, the right hand side falls 
away. Nonlocal bulk effects provide driving terms that are like anisotropic stress terms 
in general relativity In the latter case however, the evolution of anisotropic stress is 
determined by the Boltzmann equation or other intrinsic physics. 

What we can conclude from Eq. ( |73D is that the local bulk effects enhance tensor pertur- 
bations for non-inflationary expansion, and suppress them during inflation, since inflation 
imphes p + 3p < (at all energy scales) by Eq. (p!7|). The nature of the nonlocal bulk 
effects carried by P^i, requires knowledge of the off-brane dynamics of Sab, which we have 
not considered. 



E. Cosmological vector perturbations 

The linearized vorticity propagation equation (|50|) does not carry any bulk effects, and 
vorticity decays as in standard general relativity. However, the gravito-magnetic divergence 
equation ( ^8]) shows that the local and nonlocal energy density and the nonlocal energy flux 
provide additional sources for the brane gravito-magnetic field: 

D^H^, = k\p+ p)uj^ + - [{«V(P + P)+ 8W} - 3curl . 

In standard general relativity, it is necessary to increase the angular momentum density 
K^{p-\-p)u)^ in order to increase the gravito-magnetic field, but bulk effects allow an increased 
gravito-magnetic field without this. In particular, unlike in general relativity, it is possible to 
source vector perturbations even when the vorticity vanishes, since curl may be nonzero. 
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VI. CONCLUSION 



By adopting a covariant approach based on physical and geometrical quantities that 
are in principle measurable by brane-observers, we have given a comprehensive analysis of 
intrinsic cosmological dynamics in Randall- Sundrum- type brane-worlds. This has allowed 
us to carefully delineate what can and can not be predicted by brane observers without 
additional information from the unobservable bulk. 

Our main result is probably that scalar perturbations on super-Hubble scales can be 
evaluated intrinsically on the brane. This is a generalization of the result given in i.e., 
the evaluation of adiabatic scalar perturbations on super-Hubble scales during inflation on 
the brane, done under the assumption that may be neglected. 

We showed that tensor perturbations can not be evaluated intrinsically on any scales. 
This is not surprising, since the ability of gravitational waves to penetrate the 5th dimension 
inevitably introduces unpredictability from the standpoint of observers confined to the brane. 
Vector perturbations on the brane can be generated even in the absence of vorticity, via the 
curl of the nonlocal energy flux. 

Our nonperturbative results include: 

(i) Calculating the gravitational and non-gravitational acceleration felt by brane observers, 
allowing us to provide covariant characterizations of gravity and matter localization, and 
showing in particular that the non-gravitational off-brane acceleration of fundamental world- 
lines is towards the brane during inflation. 

(ii) Showing how bulk effects can disrupt the relation between isotropy of the CMB and 
spatial isotropy and homogeneity on the brane. 

(iii) Showing how bulk effects modify the dynamics of gravitational collapse and singularity 
formation. 

We have derived the exact nonlinear equations governing cosmological dynamics on the 
brane, and the corresponding covariant and gauge-invariant linearized equations. Further 
implications of these equations could usefully be pursued. In particular, an important topic 
for further research is the calculation of scalar perturbations on very large scales, and the 
investigation of limits imposed by observations. 

However, the major further step required, and not undertaken here, is to complete the 
picture by investigating the dynamical equations of the gravitational field off the brane. A 
starting point is provided by the general equations given in [^p!0|, which determine 

^n^AB , ^iiBaBC , ^iiRaBCD , 

where Rabcd is the 4-dimensional Riemann tensor, and 

SaBC = QA^ QB^CoECFn^ ■ 

However, it may turn out to be more useful to develop an alternative decomposition of the 
bulk Weyl tensor, along a timelike direction rather than a spatial direction n-^. Intuitively, 
this may provide a more direct and transparent route to the evolution equation for Vfj,^-, whose 
absence leads to unpredictability on the brane. Such a timelike decomposition requires a 
generalization to higher dimensions of the 4-dimensional decomposition of the Weyl tensor 
ri| . The generalization has been given in [|H . 
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The complete and closed system of dynamical equations would allow us to develop more 
systematic and probing tests of the Randall-Sundrum-type models against observational 
constraints. Despite the appealing geometric and particle-physics properties of such mod- 
els, it is their confrontation with cosmological observational tests that will provide a more 
decisive arbiter as to whether they are viable generalizations of Einstein's theory. 



Note: Since this work was completed, a number of papers has appeared, setting up the 
5-dimensional formalism (metric-based) of bulk perturbations p2[, and gravitational waves 



produced during inflation on the brane have also been studied 
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APPENDIX A: GENERAL PROPAGATION AND CONSTRAINT EQUATIONS 



For a general, imperfect energy- momentum tensor, as in Eq. (^, the propagation and 
constraint equations (p8D-(^) are generalized to: 

Propagation: 



1 fK 



12 \k 



^{p) + f + |curl A^ - a^^u" = 



(Al) 
(A2) 



1 



24 V/t 



(A3) 



^(^i.) + QE^^ - curliJ^j, + {p + p)(y^^ + + D(^g^) + g/t Gtt^^ + k Ai^^q^) 

1 



72 V/t 



-k'' {3 (2p^ + 2pp - TTaf^ir'"^ - 2g„g") (x^^ + 3(p + ]3)7r^^ + 3(p + 3p)7r(^,,) 

- 67r„(^7r^)" - 6g(^g^) + |D(^, (vr^)^?" - 4pg^)) + 9 ([p + 3p]n^,^ - tt^^^tt^)" - 

- 3v4(^ (4pg,.) - 7r,,)ag") + 3(p + 3p)a°' {f,TT^)a - 3o-"(^ (jc^ (u)rra)f3 + q^)qa) 

+ 3{p + 3p)uj°'eai3{,,Truf - 3t^o£:"'^(/, (tt^^^)??^)^ + } - 48W(t^i. - 36P(^,,) 



- 36D(,,Q,) - 126^^, - 72A(^Q,) - 36a"(^P,)„ - 36^"£„^(^n 



/3 



(A4) 



+ QHnu + CUrlE.i, - }:K^Cm\7Tau " 3(Ja(uHu)" + U}"" ScBtuHu)^ 



+ 2A°'ea|3{^lEu)'^ + ^K^uji^^qy) — ^fi~a"'(^^e^)apq'' 



1 ,.2„a 



1 



48 V/t 



{-2curl ((p + 3p)7r^^ - 71^,(^71^,)" - - 12pu;(^g,,) + 4cj(^7rj,)„g° 



+ (4pg^ - vr^^g^)} + 24curiP^, - 72a;^^Q,) + 24a"(^£,)„^Q'= 



(A5) 



Constraint: 



D^^^ - A'^uj^ = , 



(A6) 



1 /K 



DX,-curla;^-|D^e + «:V + 2[^,^]M = -7T7 - t^^^pq^,- ^^^^q") +2AQ^ , (A7) 



24 \k 



curl cr^^ + D(^u;,,) - ff^^ + 2A(^^Uy) = 



(AS) 
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1 / \ ^ 

+ |pD^p - 47r°^D^7r,^ + a^, {Apq^ - 7r^"g,) + 3£^„^a;°7r^^g^ - 4p[a;, g]^} 

+16D^W - leeS^ - 240^^^^, + 24a^,Q^ - 72[a;, Q]J , (A9) 

1 / ~ \ ^ 

= ^ (^^ j {curl (TT^.g'^ - 4pg^) + 4 (2p2 + 2pp - vr^^vr"^ - 2g,g") 

- 2£^„^(T% (tt.^^tt^)'^ + g<^g^>) + 2 ((p + 3p)7r,,, - tt.^^tt,)" - g^^g,)) 

+ 2(p + 3p) [or, vr]^ - 24curl + 64Wcu,, - 24[a, P]^ - 24P^,cu'^] . (AlO) 



The 4-dimensional general relativistic results are regained by setting all the right hand 
sides of these equations to 0. 



APPENDIX B: DIFFERENTIAL IDENTITIES 

On a flat Friedmann background, the following covariant linearized identities hold pO| : 



D^/=(D^/)-+/7D^/-M^, (Bl) 

D2(D^/) = D^(DV) + 2/^^, (B2) 

(D^/)- = - 2iJDV + fWA^ , (B3) 

curlD^/ = -2/cu^, (B4) 

curlD(^D,)/ = 0, (B5) 

(D^K)' = D^K-ifD^K, (B6) 

D[^D,]K. = = D[^D,]1^,^, (B7) 

D'^curl V;, = , (B8) 

D'^D(^K) = |DV^ + ID^(D'^K) , (B9) 

curl D(^K> = |D(^curl K) , (BIO) 

curl curl = -DV^ + D^(D^K) , (Bll) 

(D^iy„^)- = D^iy„^ - HB^W^p , (B12) 

D'^curl W^, = icurl (D"!^^,) , (B13) 

curl curl W^, = -"D^W^, + |D^^D"M/,)„ , (B14) 



where = V(^) and W^i, = Wi^^y) vanish in the background. 
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